Search for Spin-Dependent Short-Range Force Using Optically Polarized ^He Gas 
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We propose a new method to detect short-range P- and T- violating interactions between nucleons, 
based on measuring the precession frequency shift of polarized "^He nuclei in the presence of an 
unpolarized mass. To maximize the sensitivity, a high-pressure '^He cell with thin glass windows 
(250 /^m) is used to minimize the distance between the mass and '^He. The magnetic field fluctuation 
is suppressed by using the ^He gas in a different region of the cell as a magnetometer. Systematic 
uncertainties from the magnetic properties of the mass are suppressed by flipping both the magnetic 
field and spin directions. Without any magnetic shielding, our result has already reached the 
sensitivity of the current best limit. With improvement in uniformity and stability of the field, we 
can further improve the sensitivity by two orders of magnitude over the force range from 10"* — 10"'^ 
m. 

PACS numbers: 13.88.-|-e, 13.75.Cs, 14.20.Dh, 14.70.Pw 
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The possible existence of new forces with weak cou- 
plings and macroscopic ranges have been proposed by 
several authors [l|-y] . A P- and T- violating macroscopic 
force with an interaction range from cm to /xm first pro- 
posed in b\ has the form 
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where gs and gp are the scalar and pseudoscalar cou- 
pling constants, fi is the Plank's constant, a is the spin 
of the polarized nucleon, f = r/r is the unit vector from 
the unpolarized nucleon to the polarized nucleon, 7ti„ is 
the nucleon mass, and A is the range of the force. This 
short-range force is mediated by exchanging an axion- 
like particle between unpolarized nucleons and polarized 
nucleons. A similar interaction may also exist between 
nucleon and electron. Many experimental efforts have 
been devoted to search for this interaction between either 
nucleons or nucleon and electron, and various techniques 
have been used, such as sensitive torsion pendula [y, |7J, 
clock comparisons between two different polarized species 
[8|-[l0j , and measurements of neutron bound states on a 
flat surface in the gravitational field [ll|. Very recently, 
measurements of the longitudinal relaxation rate Ti and 
transverse relaxation rate r2 of polarized ^He gas were 
used to search for this short-range interaction between 
nucleons 12h15|. As the relaxation time of polarized 
•^He can be as long as tens of hours, any new interac- 
tion with the polarized '^He nuclei can lead to a visible 
change in the relaxation time. These measurements pro- 
vide to our knowledge the most stringent direct labora- 
tory constraint on the coupling constant product gsgp for 
a monopole-dipole interaction between nucleons of the 
form in Eq. ^ over distances from 10~^ to 10~^ m 'l5j. 
Note that this limit is still more than 9 orders of magni- 
tude larger than the standard Axion coupling originally 
proposed to solve the strong CP problem [8|, ll6| . 



In this work, we present a new method to search for 
the spin-dependent macroscopic force between nucleons 
by measuring the frequency difference of optically po- 
larized ^He gas with and without a nearby unpolarized 
mass. The frequency difference due to the magnetic field 
gradient is a first order effect, as such it is more sensitive 
than the relaxation measurement because the gradient- 



induced relaxation is a second order effect [171 , Il8l | . We 
also performed a pilot experiment to demonstrate how 
this method works. With a modest stability of the mag- 
netic field, the sensitivity of this experiment already 
reaches the current best laboratory limit on gsgp- With 
improved stability of the magnetic field, the proposed 
method could be used to improve the current best limit 
by two to three orders of magnitude in the force range 
from 10""' to 10~^ m. 

The spin-dependent short-range interaction changes 
the precession frequency of the polarized nuclei through 
the interaction ct • f in Eq. ([T]), which is similar to the 
well-known fi-B interaction of a magnetic dipole moment 
/t in an external magnetic field B. Consider a cylindri- 
cal cell containing polarized '^He gas with its polarization 
pointing in the z direction along the axis of the cylinder 
and a block of unpolarized mass is placed next to the end 
of the cell. The short-range interaction on each '^He atom 
inside the cell can be obtained by integrating Eq. ([T]) over 
the unpolarized source mass. In the limit case in which 
the transverse dimensions (x, y direction) of the mass are 
much larger than the force range, the mass can be ap- 
proximated as an infinite plane source with its normal 
pointing in the z direction. In this limit the frequency 
shift from the planar mass block is [14] 
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where z is the distance measured from the surface of the 
mass block to the polarized "^He, N is the nucleon number 



density of the mass, and do is the thickness of the mass 
block. If the external magnetic field is uniform with a 
value of Bo, then the spin precession frequency in the 
presence of the mass has a spatial dependence which can 
be written as 
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) . As the precession signal 



received by the pickup coil is a weighted sum from all the 
■^He inside the cell, the signal induced in the pickup coil 

is 



cos 



(wot + ^e"^/^t)B(z)dz, 
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where d is the window thickness of the cell, B{z) is the 
field profile of the pickup coil along the cell axis, and the 
reciprocity theorem is applied here to compute the signal 
induced in the pickup coil [19]. When A is zero, Eq. (|4]) is 
a pure sinusoidal function with a well-defined frequency. 
When A is nonzero Eq. ^ shifts the mean frequency of 
the signal. The mean oscillation frequency determined 
from Nc observed periods during a time T is f = Nc/T. 
In the presence of the interaction for the same number 
of periods the time duration changes to T' and the new 
frequency is /' = Nc/T'. The frequency difference is ^ 






(5) 



where AT ^ T' -T. Eq. ^ establishes a relation be- 
tween A/ and AT. AT can be calculated for different 
strengths gsQp and ranges A of the spin-dependent in- 
teraction by numerically integrating Eq. Q. For any 
given A, the parameter A in Eq. Q is tuned in such a 
way that the calculated AT matches the experimentally 
determined frequency shift A/. Hence, constraints on 
gsQp with different values of A are established given the 
sensitivity of the experiment. 

In this pilot experiment we used a 7 amg ^ high pres- 
sure '^He cell originally constructed as a ^He gas tar- 
get for two- and three-body photo-disintegration exper- 
iments [20|. The cell had two chambers, a spherically- 
shaped spin-exchange optical pumping chamber and a 
40 cm long cylindrical target chamber connected by a 
glass tube (see Fig. [T]). The target chamber had two 
thin glass windows on its ends. The thickness of the win- 
dow was about 250 /xm. A Macor machinable ceramic 
mass block was used as the unpolarized mass. It was 
repeatedly brought into contact with and moved away 



^ When S/N becomes the limiting factor of the measurement (not 
the case in the present work) , phase difference or frequency spec- 
trum may provide a better measurement of the frequency shift 
than the peak counting method. 

■^ 1 amg is the number density of 1 atm gas molecules at °C. 



from the cell window by a G-10 rod connected to a step- 
per motor. The stepper motor moved the mass to a final 
position with better than 10 micron repeatability, more 
than an order of magnitude smaller than the cell window 
thickness. Two identical pickup coils were mounted right 
below each window. Pickup coil A was used to measure 
the frequency of the polarized ^He nuclei infiuenced by 
the mass; pickup coil B was used to monitor the magnetic 
holding field. The 40 cm long '^He cell was positioned at 
the center of a Helmholtz coil pair. Due to fairly large 
gradients at the end of the cylindrical chamber, two iden- 
tical two-axis gradient coil pairs (dashed line in Fig. [1]) 
were added at both ends to compensate the gradients 
from the Helmholtz coils and background fields. The 
measured transverse relaxation time T2 of the '^He signal 
is greater than 1 s with this arrangement. 
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FIG. 1: The diagram of the test experiment apparatus (not 
to scale). The cylindrical cell axis is the z direction. The 
cell contains 7 amg He gas and is optically pumped in the 
pumping chamber to about 40% polarization. The coils in the 
dashed lines are gradient coils to actively compensate the gra- 
dients from the Helmholtz coil and other background fields. 

Free induction decay (FID) at 24 kHz was performed to 
measure the '^He precession signal. The RF pulse with 
small tipping angle was applied to make the polariza- 
tion loss negligible. The precession signal is digitized and 
recorded by the computer. In order to determine the fre- 
quency unambiguously, the acquisition time stopped at 
the instant when the signal-to-noise ratio is either below 
10 or at 0.2 s, whichever comes first. The frequency is 
computed by counting the periods during the acquisition 
time. In a 7 amg ^He gas cell, the diffusion constant is 
about 0.27 cm^/s [21|, and it is known that the effective 



diffusion rate is lowered at the cell boundary [22[ . There- 
fore effects from the diffusion of the '^He can be ignored 
during the FID measurement and it is valid to use Eq. 
(|H) to analyze the experimental data. Each measurement 
cycle contains two measurements: first with the mass in 
contact with the window (in position) and next with the 
mass moved far away from the window (out position). 



Simultaneous measurements at pickup coil B were per- 
formed to monitor the magnetic field fluctuations. The 
peak to peak variation of the field is about 5 x 10~'^%. Re- 
moving the field fluctuations measured by coil B reduced 
the peak to peak variation of the field to 4 x 10""'%. After 
the field correction, the frequency difference between the 
"in" and "out" positions is calculated as A/ = fin — font- 
The magnetic susceptibility of the mass can change the 
field at the location of the ^He through its effect on the 
holding field and therefore lead to a systematic effect. Al- 
though the magnetic susceptibility of Macor ceramic is 
known to be small enough not to cause a systematic error 
in this test experiment, the real material can in principle 
contain paramagnetic or even ferromagnetic impurities. 
Paramagnetic impurities would increase the local field 
strength (and therefore the '^He precession frequency) in- 
dependent of field direction. The spin-dependent inter- 
action can increase or decrease the precession frequency 
depending on the magnetic holding field direction. We 
therefore can isolate a possible spin-dependent interac- 
tion from paramagnetic effects by flipping the magnetic 
field. However, a frozen-in field from possible ferromag- 
netic impurities has the same magnetic field dependence 
as the spin-dependent interaction, which makes it dif- 
ficult to separate them apart. The most likely ferro- 
magnetic contamination of the mass block comes from 
machining process, during which ferromagnetic tools are 
usually used to cut the material. To minimize this ef- 
fect, cutting tools with diamond tips are used to ensure 
that there is no physical contact between the metallic 
part of the tools and the surface of the block. Ad- 
ditionally, we also flip spin direction in order to can- 
cel any spin-dependent systematic effect in the system 
if there is any. Therefore, we took data in four differ- 
ent configurations for B field and spin directions, A/++, 
A/_+, A/__, and A/+_, representing B/S: ++, B/S: 
— 1-, B/S: -, and B/S: -|— , respectively. If A/s represents 
the field-dependent frequency shift, A/5 represents spin- 
dependent frequency shift, and A/o represents frequency 
shift without B or S dependence, then 



A/++ = 
A/-+ = 
A/- = 
A/+- = 



-A/b + A/5 + A/o 
-A/b + A/5 + A/o 
-A/b - A/5 + A/o 
-A/b - A/5 + A/o. 



(6) 

(7) 
(8) 
(9) 



In this notation, the short-range force induced A/b can 
be expressed as 

A/b = i(A/++ + A/+_ - A/__ - A/_+). (10) 
The uncertainty of A/b is given by 

4 

where (7+-|_, a y, a , and a^ are the uncertainties 

of A/++, A/_-|_, A/ , and A/^ , respectively. The 
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noise in the measurement is mainly due to the magnetic 
field fluctuation, which limits the uncertainties of the fre- 
quency measurements. 

We performed 100 measurement cycles for each of the 
four configurations to determine the average frequency 
difference between the two mass positions. The frequency 
differences from the spin- dependent force for these runs 
are shown in Fig. [5] The average frequency difference 
of 100 measurements was A/b = -0.003 ± 0.005 Hz, 
consistent with zero. 
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FIG. 2: (Color online) The frequency difference correlated 
with the position of the ceramic mass block. The error bars 
show the standard deviation of the magnetic holding field 
after correction by coil B. 



Eq. ([T0|) is then used to calculate the frequency differ- 
ence from the spin-dependent short-range force to place 
an upper limit on gsgp as a function of A. For / = 24 kHz 
and A/ = 0.005 Hz, Eq. © yields AT = -2.1 x lO^^T s. 
With a fixed frequency shift, AT should increase linearly 
with respect to the acquisition time T. Theoretically, AT 
is obtained by comparing Eq. (j4|) with a sinusoidal func- 
tion in absence of the force. With a real pick-coil profile, 
Eq. ^ is obtained by numerical integration, using the 
actual geometry of the experiment. The resultant AT as 
a function of T is shown in Fig. [3] (the upper black curve). 
Surprisingly, AT increases linearly only for a short pe- 
riod of time. As time elapses, AT oscillates around a 
constant value, indicating that the frequency shift due to 
this exponential type of force is not fixed and diminishes 
at large T. This striking behavior suggests that one will 
not gain more information from longer measurement time 
though the frequency resolution is improved by doing so. 

An closed-form solution of Eq. (j4]) can be obtained if 
one approximates the real profile oi B(z) by a rectangular 
function with a cut-off position w mimicing the width of 
the profile. In this case, the upper limit of the integral is 
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FIG. 3: (Color online) The upper curve shows the |Ar| as a 
function of T, using the real field profile of the pickup coil. 
The lower curve shows the same curve but with the pickup 
coil profile approximated by a rectangular function. The inset 
of the figure shows the linear behavior of | AT] at small T. 



replaced by w, and the integration yields 

S{t) = \[cos{uont){Ci{Ae-it) - Gi{Ae-'^t)) 
- sin(wot)(Si(Ae~^t) - Si(yle"^t))], 



where Ci(x) is Cosine Integral and defined as Ci(x) = 
— /^ ^^dt; Si(a;) is Sine Integral and defined as Si(a;) — 
J^ ^dt [Hi. One can use Eq. dH]) to calculate the AT 
in Eq. ([5|). To the first order, AT can be expressed as 



AT: 



SijAe-T^T) -Si( Ae-fT) 

Ci{Ae-iT) - Ci{Ae-fT)' 



(13) 



Using the value oi d = 250 /xm and w = 5.825 cm, Eq. 
p3)) is plotted in Fig. |3] (the lower red curve). It is sur- 
prising to find that AT does not increase linearly with 
respect to T all the way up. At certain point, AT be- 
comes more or less a constant, which means the phase 
diff'erence due to the short-range force stop accumulat- 
ing after certain time. A longer measurement time will 
not increase the phase difference due to the short-range 
force. This behavior is the result of the inhomogeneous 
broadening due to the varying distance between the '^He 
spin and the mass. 

This simple approximation yields a satisfactory result 
in terms of the time at which the linear relationship 
breaks down (The exact solution is plotted as a black 
curve in Fig. [3|). In our experiment, the maximum mea- 
surement time is less than 0.2 s to guarantee that the 
FID measurement is in the linear region, so that the fre- 
quency comparison is valid even the measurement time 
for each configuration is slightly different. 

By choosing different values of A, the constraints on the 
coupling constants gsffp are found and plotted as a solid 
line in Fig. S) The 250 /im window thickness of the dou- 
ble chamber glass cell allows us to constrain interactions 



ranges A down to 10~^ m. The dominant source of the 
uncertainty in our experiment came from magnetic field 
fluctuations. The gradient compensation coils, needed to 
achieve an uniform field in our apparatus, added uncor- 
related magnetic field noise to the holding field as the 
gradient coils were powered by independent power sup- 
plies. 

There are many avenues for the improvement of the 
measurement sensitivity using this technique. One can 
use a dedicated ^He cell with a shorter length and add 
magnetic shielding instead of gradient coils to improve 
the field uniformity and stability. In [2J], the authors 



conducted a precision frequency measurement using po- 
larized gases in an apparatus with three layers of cylin- 
drical fj, metal shielding for field uniformity and a co- 
magnetometer technique to reduce the magnetic field 
noise by 3 to 4 orders of magnitude. They achieved a pre- 
cision of 10~^ Hz out of 1000 Hz, two orders of magnitude 
better than our pilot experiment. In Fig. |4]we show the 
limits on gsgp (dotted line) which could be achieved with 
this technique using the sensitivity demonstrated in [2J] 
assuming all other geometric parameters (window thick- 
ness, pickup coil size, etc.) are unchanged. As shown in 
FiglH this projected sensitivity would represent a signif- 
(12) leant improvement in the A region of 10~^ to 10~^ m. 



compared with the best existing laboratory limit derived 
by the authors of [l5|, shown as the dotted-dashed line, 
based on the results of a T2 measurement from [25| . 

The sensitivity of the experiment below A = 10~^ m 
is clearly limited as the thickness of the window becomes 
much larger than the force range. The 7 amg ^He cell 
has an internal pressure of more than 10 atm at 200 °C 
and the 250 /xm window thickness is needed for strength. 
However, the high pressure is not necessary for this ex- 
periment. A 1 amg '^He cell with reduced window thick- 
ness could be used to improve the sensitivity of the mea- 
surement, especially in the A < 10^"* m range. Another 
order of magnitude improvement on the sensitivity could 
be achieved if denser material was used as the unpolar- 
ized mass, such as pure copper or tungsten. We conclude 
that our proposed method shows a promising sensitivity, 
with at least one to two orders of magnitude improve- 
ment over the current best limit possible in a dedicated 
experiment with better magnetic field stability. A even 
higher sensitivity could be achieved if a thinner wall for 
the "^He cell and denser material is used. 
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FIG. 4: (Color online) Constraints on the coupling constant 
product gsQp of the spin-dependent force as a function of the 
range A and the equivalent mass of the axion-like particle 
mediating the short-range interaction. The dashed line is the 
result from [3|, the dash-dotted line is the re-analysis of the 
T2 measurements of [23| by [I3I1 the solid line is the analysis 
of our present experiment, and the dotted line is a projected 
sensitivity achievable using our method based on the stability 
of the magnetic field demonstrated in [241 ]. The dark gray is 
the excluded region and the light gray is the region that could 
be excluded with the improved field stability. 
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